Two-point function of strangeness-carrying vector-currents 
in two-loop Chiral Perturbation Theory 



Stephan Diirr*"'^ and Joachim Kambor^'^ 
"■University of Washington, Physics Department, Seattle, WA 98195-1560, U.S.A. 
^Paul Scherrer Institute, Particle Theory Group, 5232 Villigen, Switzerland 
University of Zurich, Institute for Theoretical Physics, 8057 Zurich, Switzerland 



Abstract 



We calculate the correlator between two external vector-currents having the 
quantum-numbers of a charged kaon. We give the renormalized expression to 
two loops in standard chiral perturbation theory in the isospin limit, which, 
as a physical result, is finite and scale-independent. Applications include a 
low energy theorem, valid at two loop order, of a flavor breaking combination 
of vector current correlators as well as a determination of the phenomeno- 
logically relevant finite 0(p^)-counterterm combination Qy by means of in- 
verse moment finite energy sum rules. This determination is less sensitive to 
isospin-breaking effects than previous attempts. 



Preprint numbers: UW-PT/99-14, ZU-TH/15-98, [hep-ph/990753g 
PACS numbers: 12.39.Fe, 11.55.Hx 



*stephan.duerr@psi.ch 
^kambor@physik.unizh.ch 







I. INTRODUCTION 



The appropriate tool to analyze hadronic processes at very low energies is chiral pertur- 
bation theory (XPT), the low energy effective theory of QCD. The general structure of XPT 
is well understood and it has been successfully applied at the one- and two- loop level 

to a vast body of low-energy data. Recently, the non- anomalous counterterm Lagrangian 
at Oijf') has been constructed in full generality for chiral SU{N) thereby closing a gap 
left behind in an earlier attempt The divergent part of the generating functional at the 
two- loop level has also been calculated in closed form and is now available [0. 

Despite this advanced state on the theoretical side, phenomenological applications at 
the two- loop level are relatively scarce []8|-|10|]. Two reasons may be responsible for this 
shortcoming: being an effective theory, XPT introduces new, unknown coupling constants 
at each order of the perturbation series. The number of such unknown parameters is rather 
large at 0{p^) Consequently, it is a priori not clear whether at this order independent 
physical processes can be related in a parameter free manner - at one-loop one of the 
main virtues of XPT. The second, more technical reason, is the large computational effort 
demanded to perform complete two-loop calculations, comprising the finite contributions to 
loop-integrals and counterterms, in particular in the 3-flavor case. 

In this paper we calculate the two-point function of strangeness carrying vector currents 
to 0{p^) in standard XPT. Although completely analogous to the calculation of the cor- 
relator of flavor diagonal currents considered in 0, the non-equal masses propagating in 
the loops make the present investigation technically more demanding. We use a method 
where one-loop subgraphs are renormalized before the second loop-integration is performed, 
which helps simplifying intermediate expressions. Combining our result with the isospin and 
hypercharge component of the vector correlator obtained in 0, we can form a "low energy 
theorem", i.e. a parameter free relation between physical quantities, valid at 0{p^). This 
shows that, despite the many coupling constants present at this order, there are still specific 
flavor breaking combinations of observables which can be predicted in terms of meson masses 
and decay constants only. 

A more phenomenological motivation of the present calculation derives from experimental 
efforts to determine the spectral functions of strangeness carrying currents. The ALEPH 
collaboration has recently published their analysis of r-decays into hadronic final states 
with strangeness This concatenation of data allows to connect, via inverse moment 



finite energy sum rules (IMFESR), the vector current correlator at low energies calculated 
in this paper to the asymptotic behavior of the correlator obtained by the operator product 
expansion |jl2|. As a result, we are able to fix a coupling constant of the 0{p^) chiral 
Lagrangian, Qy, from the difference of isovector and strange component of the vector spectral 
functions. The constant Qv has been estimated in Ref. [0 from the difference of isovector 
and hypercharge components of the vector spectral functions. However, as pointed out 
11^ , that determination is sensitive to isospin violating corrections affecting mainly the 



m 



hypercharge current. The strange current on the other hand is not sensitive to isospin 
violation, and the determination of Qv presented in this paper is therefore more reliable. 

The article is organized as follows. Section II contains an outline of the general strategy 
employed in the calculation and a derivation of the strange formfactor which plays a central 
role in it. In Section III the two-point correlator {T{VsV^}) between two strangeness carrying 
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external vector-currents is derived. The physical content of this result is discussed in Section 
IV where we also point out several consistency checks. Combining strange, isospin and 
hypercharge component of the vector two-point function yields a new low energy theorem 
valid at 0{p^). In section V we present a determination of the low energy constant Qy 
by employing an inverse moment finite energy sum rule for the difference of isovector and 
strange component of the vector current spectral functions. Finally, section VI contains our 
conclusions. Some technical details are collected in the Appendix. 



II. FORM FACTOR CALCULATION AND DESCRIPTION OF THE METHOD 

For explicit calculations we shall use the Lagrangian of standard XPT 

c = c^^^ + + c^"^ (1) 

where denotes the Lagrangian of order p\ The first two terms in this expansion can 
be found in the work by Gasser and Leutwyler 0,^] which also discusses how the meson 
fields couple to external sources. The last term, i.e. the one of order p^, has recently been 
constructed for the non-anomalous sector in 0, thereby reducing the older form given in 
to a minimal set. We shall employ the SU(3) version of the Lagrangian given in [^]. 

From (P, the correlation function between two external vector currents with the quantum 
numbers of a kaon could be computed following the usual procedure: Draw all Feynman 
diagrams contributing at 0{p'^) and 0{p^) in standard XPT [^|[ (see Fig. |] and Fig. |[), 
evaluate them and have the divergences in the loop- and counterterm-contributions canceled 
to get the finite physical result valid in 4 dimensions. Naively following this recipe would 
amount to compute the two-loop integrals in terms of unrenormalized masses and couplings. 
As it turns out, the two- loop contributions are significantly shorter if they are expressed in 
terms of one-loop renormalized expressions. 

A well known method to do this is simply combining diagrams (e) and (f ) in Fig. |^ into a 
new single diagram, equivalent in topology to the one-loop unitary diagram, employing mass- 
and wavefunction-renormalization correct to O(p^). Likewise, diagrams (b) and (c) can be 
replaced by a single new diagram, topologically equivalent to the one-loop tadpole-diagram, 
without altering the expression at 0{p^), if one-loop masses are used. In the remaining six 
diagrams in Fig. ^ the K, vr, //-masses may be taken as the one-loop renormalized masses as 
well, since the corresponding change shows up at 0{p^) only. 

An alternative method uses the strategy of renormalizing intermediate results in an even 
more direct, physical way. It is known as the "master method" and was introduced in 
the context of XPT in ref. [[l^]. The idea is to start from the one-loop expression for the 



formfactor of the strangeness-carrying vector-current with any of the possible intermediate 
particle-combinations of the V^V^-correlator in the final state. The correlator follows by 
contracting this formfactor with its hermitian conjugate (modulo a modification necessary 
to avoid double-counting, c.f. Sec. Ill), integrating over internal momenta, adding tadpole- 
and counterterm-contributions and completing the remaining renormalization procedure. 

We have calculated {T{VsVj}) following either one of these two methods and have found 
the results to agree. For this reason we restrict ourselves to presenting the calculation in 
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the second, more physical framework. Below we summarize the results for the strangeness- 
carrying formfactor, while section III contains the derivation of the correlator. 

It is clear that the technical key quantity needed in the second strategy is the off-shell 
form-factor S^[PiP2\{pi,p2)i i.e. the amplitude for an external vector-current with momen- 
tum = {pi +P2)fi and the quantum numbers of a to decay into any of the combinations 
P1P2 = {K~^n^ , K'^T], K^n'^) with arbitrary momenta pi and p2, respectively. Note that S*^ 
is not a physical quantity, i.e. the result depends on the choice of interpolating fields for 
the Goldstone bosons. However, the correlator we are aiming at involves an integration 
over these fields and is independent of this choice. For convenience, the formfactor shall be 
separated into a symmetric and an antisymmetric piece w.r.t. the outgoing momenta 

S^[Pi P2] = S4P, P2] {pi + P2)^ + 5- [Pi P2] {Pi - P2)m (2) 

and the convention is that pi denotes the momentum of the strange meson (Pi = K~^, K^), 
whereas p2 is the momentum of the non-strange meson (P2 = vr^, 77, vr"*"). Fig. ^ displays 
diagrammatically the relevant contributions. Diagram (a) represents the leading 0{p'^)- 
contribution, whereas diagrams (b, c) as well as (d, e, f) are of order 0{p'^). The latter three 
are frequently referred to as the "unitary"-, "tadpole"- and "counterterm" -contributions. 
The sum of these six diagrams corresponds to a single formf actor-diagram in the effective 
theory [in the field-theoretic sense, i.e. in the theory where all quantum corrections have 
been absorbed into the propagators and vertices of the effective action and Greens functions 
are constructed from tree-level diagrams only] where the two outgoing mesons carry the full 
wavefunction-renormalization factors Zp-^, Zp^, with Zp = 1 — Ap and 

A. = - ^^(O) + 4((2< + M'^L, + MIL,) (3) 

-f^o 3 3 i^o 

= jr2{-\A{Ml) - \a{MI) - \a{M'4) + ^((2Mi + Ml)L, + MIL,) (4) 

A, = '{-A{Ml)) + ^{{2Ml + Ml)L, + MJLs) . (5) 

However, in the light of the application to follow, it is extremely convenient to have the 
external legs carry a factor a/Zp, rather than Zp or 1. Since Zp = 1 — Ap implies y/Zp = 
1 — ^Ap + 0{p'^) and since we are interested in to 0{p'^), this "-\/Zp-convention" amounts 
to defining S'^[PiP2] as a weighted sum: Diagrams (a, d, e, f) are included with their full 
weight, but diagrams (b, c) with a factor | only. Note that renormalizability of the on-shell 
formfactor is not affected by such a reweighting as long as diagrams (b) and (c) receive the 
same factor and the group (d, e, f) also receives a common factor. In practice, on-shell 
renormalizability is implemented as follows: First replace pf — > M"^, p\ Mp^ (P2 = 
7T~^ , 71^ , rj) , where Mp denotes the one-loop mass of particle P. Next, expand all integrals 
in — 4), or equivalently in A, as is explained in the Appendix, and simultaneously the 
renormalization constants according to 

L, ^ /-^ ■ (Lf^A + + 0(A"')) . (6) 

Finally, check that there is indeed a set of numerical constants {-^vj^''} for which all contri- 
butions in A cancel exactly. These "divergent" parts of the one-loop counterterms agree of 
course with their standard values l\^^ = [@,|^. 
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In summary, the strange off-shell formfactor takes the form (0) with 
-2iLf\Ml-Ml)+iLf{p\-pl) 



+ l[Mi + MI + 5{pl +pI- 3q^)] B{q\ Ml) 



1 M"^ — A/f^ 



+ 



3 vl-vl 
16 g2 

-1 



[(Ti? - T«)(g^ Mi, M^^) + {T^2 - T^2){q\ M^, M,^)] } (7) 



2iLf\Ml-M^^)+iLf\p\-pi) 
l[Mi - 7M^^ - 3{pl +pI- 3g^)] ^^^^^ 

[19Mi - + 3(p? + - 3g2)] ^ ^ " M^, M^^ 



■-[19Mi - 5M,2 + 3(p? + _ 3g2)] " , ...^ , ...^^ 



+^ ^[(T^ - T«)(g^ Mi, M.^) + (T^ - T«)(g^ Mi, M,^)] } (8) 
5+[ir°7r+] = v^SjK+7r°] (9) 



3 

16 



[Afi,(M2) + 2Afi,(Mi) + AUMl)] - i4°^g' 



(10) 

(11) 
(12) 



where the particle masses Mp refer to the one-loop renormalized masses and the integrals 
are all finite and explicitly given in the Appendix. 

It is worth mentioning that this result reduces to the expression known in the literature, 
if we restrict ourselves to on-shell outgoing momenta. Employing the notation of 



{K^,p'\V2u^,s\7i\p) = {p'^+p^)f^-{t) + {p'^-p^)f^-(t) 



(13) 



and 



fP^P^{t)^f^^^-{t) + 



P1P21 



f'^'^t) 



Ml - M| 



(14) 



our result simplifies in the on-shell limit to (|T3|, |T^) with 
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= 1 + Ml, Ml) - A^M) - A^^{Ml)] - iq'L^ + {tt ^ v)} (15) 

fo^iq') = ^ + ^2^8 Ml-M^ Mfin(M^2) _ - 3Afi,(M^^)] +ALf\' (16) 

+ \ [5g^ - 2{Ml + Ml) - ^(^K- Ml)^ ^ 1-^^^,^ 
o q' 1 

fo^'iq') = 1 + {-| m^Tm^ [^fin(M,2) - 2Afi„(Mi) + Afi.(M,2)] + 4Lf (17) 

+ 1 [3g^ - 2(M^ + Ml) - ^^^~^-^' ] -5(g^ M^, Ml) 
o q' 1 

+ 1 [-9g^ + 2(Mi + 9M,^) - ^^^^ ~ ^^'^' ] i5(g^ M^, M^^) } . 
This coincides with the result given in 0. 



III. TWO-POINT FUNCTION CALCULATION 

Here we give the derivation of the technical key quantity of this paper, the correlator 
{T{VsV^}) between two external vector currents with the quantum numbers of a K^, respec- 
tively. The convention is that Vg creates a multi-particle state with the quantum numbers 
of a , or annihilates a state with the quantum numbers of a K~ , i.e. 

K" = (18) 

with Xa the standard Gell-Mann matrices and q = {u,d,s). Defining the momentum-space 
representation of the correlator by 

%.(?;., g.) = ijd'x e''^-^{0\T{Vnx) K;(0)t}|0) , (19) 

and decomposing the latter according to 

T^vAQ^^. 1^) = - q'g^^u) n'^liq') + q^^qu u'vliq'') , (20) 

the result will be given in terms of 11^^^ and 11^^^. 



A. Calculation 

For convenience we shall consider the "unitary", the "tadpole" and the "counterterm" 
contributions to illy ,, separately, i.e. the latter is decomposed as 
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mv,s = in™ + in*'^^ + in™* , (21) 

where each term on the r.h.s. shall contain both the contribution at O(p^) and at 0{p^) 
in standard chiral perturbation theory. In diagrammatic language the two-loop "unitary" 
contribution in^" is given by the sum of diagrams (e, f, g, h, i) in Fig. |^. The "tadpole" con- 
tribution in^e*^ denotes the sum of diagrams (a, b, c, d), and the "counterterm" contribution 
in^e* is just diagram (j) in Fig. ||. 

The key ingredient in the approach to the correlator {T{VsVj}) we are describing is the 
observation that the "unitary piece" inp4^p6 can be computed from the renormalized strange 
formfactor obtained in the previous chapter. The prescription is: Contract the formfactor 
S^[Pi, P2]{pi,P2) with its hermitian conjugate and sum over all possible intermediate particle 
Pi, P2 and momenta pi,P2- In practice, several points have to be paid attention to. 

First, we should mention why we defined the formfactor such that outgoing mesons carry 
only the square-root of the usual wavefunction renormalization factor, but neither Z nor 1: 
This choice allows us to proceed by splicing two of these formfactors together without having 
to worry again about wavefunction renormalization - the internal propagators generated in 
this process will automatically carry the full factor Z. The second point is that splicing the 
formfactor with its hermitian conjugate will yield the "eyeglass" -diagram (i) in Fig. |^ twice 
instead of only once. Thus one has to subtract half of the "unitary" diagram (d) in Fig. |^ 
from the formfactor before doing the splicing with the resulting object. In summary, the 
unitary piece of the correlator - which shall include the contributions both at 0{p'^) and at 
0{p^) in the chiral expansion - is constructed as 



d^i S,[Pi,P2]-S,[Pi,P2]ii) 



V+P« J (2vr)4 (£2 - Mf )((g - i)^ - Ml^ 



^ f dH ^ (2£-gV(2£-g). (gi°^' + 2^L'^-^i°^ + 0(p^))(i)^ 
(2vr)4 I (£2-Mf)((g-£)2-M|) 

{P - Ml){{q - - Mi) J" ^ ' 



In (|22D S± denotes the "modified formfactor" in which the unitary contribution is included 
with a factor |, i.e. S± = S± — \S™\ which clearly affects only the O(p^) -contributions 

in S+[Pi,P2] = 4'^[Pi,P2] + 0(/) and S_[Pi,P2] = ^i°^[Pi,P2] + 5i'^[Pi,P2] + 0(/), 
but not Si°^[Pi,P2]. Of course S± is not renormalizable, i.e. including this factor | into 
the formfactor already in the previous section is not an option. Note that this procedure 
generates an expression for inpr:^^^ where the Mp refer to the one-loop renormalized masses. 

Finally, the "unitary" piece ( P2|) must be augmented with the "tadpole" and the "coun- 
terterm" contribution. in*4'^p6 is put together from diagrams (a, b, c, d) in Fig. || plus the 
single tadpole- diagram at 0{p'^), i.e. diagram (b) in Fig. |I|. m^^^pe represents the sum of 
diagram (j) in Fig. ^and the single counterterm-diagram at O(p^), i.e. diagram (c) in Fig. |1|. 
It is important to note that both in*4'^p6 and ili^^j^^e happen to combine a 0{p'^) contri- 
bution and a 0{p^) contribution in such a way that the total "tadpole" or "counterterm" 
contribution depends on the one-loop renormalized masses only. 
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B. Renormalization 



Having added the various contributions to Spij^pS the only task left is renormalizing the 
result, i.e. rewriting it in such a way that it can be continued to c? = 4 dimensions. 

First we have to remind ourselves that the correlator is constructed from various elements, 
some of which have undergone renormalization at 0{p'^) already. In practice, this means 
that all hf^ showing up in this expression must be replaced hf^ —>■ Lf^ + \ ^ before 
the so-far untouched constants get expanded according to 

Li /-^ ■ (Lf ) A + Lf ) + L^'^X'' + 0(A (23) 

as is necessary in a two-loop calculation. Next, all integrals must be expanded in l/{d — 4), 
or equivalently in A, the latter being defined in ( [73D in the Appendix. For some more 
standard definitions in this "A-scheme" the reader is referred to [0. An explicit calculation 
in this scheme, where technical aspects are presented with great care, is the analogous 
computation of the axial correlator to 0{p^) in chiral perturbation theory |]T3|. A discussion 



of the relationship of the "A-scheme" to the MS-scheme employed in |jTO| is also found in 

The final step is checking that, after the remaining have been assigned their standard 
values Fj PJ^, there is indeed a set of numerical constants {Cj^\Cj^\ Lf^} for which all 
contributions to the correlator with a positive power of A cancel exactly. Demanding that 

— 2 

the total contribution in proportion to A is zero we find the constraints 

= 5~8{ACif + Cif)F^ (25) 

= 1 + SCifF^ (26) 

= cil^F,' (27) 

= Cg^F,' (28) 

whereas requiring the total contribution in proportion to A to be zero yields 

= 10Lf^-3iACS+Cil^)F^ (29) 

= + (30) 

= 3(Lf +lS?)-4C«Fo^ (31) 

= (4°^ + 4?)-4C«Fo^ (32) 

Constraints (22-29) prove consistent with what is known in the literature 



C. Results 



In summary, the renormalized correlator {T{VsV^}) between two external vector currents 
carrying the quantum numbers of a respectively, is found to take the form (|T9|, EH) with 
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+3iLr Mi, M^^) + T«(g^ Mi, M,^)] 

+3iLS?[Afin(M,2) + 2Afi„(Mi) + Afi„(M2)] 

-^[-5Afin(M2)2 + AA^^{Ml)A^^{Ml) + 6^fin(M2)^fi„(M2) 

+4Afin(Mi)2 - UAUMl)AUMl) + 3^fin(M^2)2] 



32g2' 

Mi, M^^) - Afin(Mi) - ^fin(M^2)]2 
■ ^^^~^-^' 0v - Pvg^ - 4Migv - 4(2Mi + M^)i?v} • 

ns ^ {-|^%^^(.^Mi,M.^) - |^^^^^(.^Mi,^)} 

+ {2iLf ; [3>lfin(M,^) - 2Afi„(Mi) - A^Ml)] 



+ ^[-bAUMlf + 4^fi,(M2)^fi,(Mi) + &AUMl)AUMl) 

+4Afin(Mi)2 - 12Afin(Mi)Afin(M^2) + ?,A^^{Mlf\ 

+^K. ■ ^^^^^icf, Mi, M^) + • ^^^^^^h icf, Mi, MJ) 

where 

KK, = ^16^fl„(M3 - 2.4j„(Mi) - 3Aa„(M,^)] - 6i(M|. - M^)L<°' 



32 q^ q^ 



32 g^ ' 



Ac;^, = -_[Afi„(M,^) - 2^fi„(Mi) + ^fin(M^^)] - 6i(Mi - M'^)Lt 



32^ ^2 — ^ gllSM^ - 3M^) + 3g ] B[q , M^, iV/J 
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In ( p3D and (|3^) the contribution at O(p^) has been separated from the one at O(p^) and 
we have used the abbreviations 



Ov = 4(4Ci? + cSVo' - (40/3)L^^) (37) 

Pv = ACil^F^ + ALt^ (38) 

Qv = ^CjS^F^ - 3{Lt^ + (39) 

Rv = ^ClSF^-iLi-'^+L[o'^) (40) 

as these combinations of (finite parts of) 0(p^)-counterterms occur quite generally 0. 
While completing this paper, we received a preprint JTOf where this result has been 



obtained independently. Up to an overall factor 2 (due to a factor -\/2 difference in the 
definition of the currents) the results agree. 



IV. DISCUSSION AND LOW ENERGY THEOREM 

In Sect. Ill we have obtained the fully renormalized finite two-loop expressions for the 
correlator of strangeness carrying vector currents. Here we shall first discuss the general 
structure of the result and point out several consistency checks which have to be satisfied 
on general grounds. Then we combine our result with the corresponding expressions for the 
isospin and hypercharge currents |^ to obtain a low energy theorem valid at 0{p^). 



A. Consistency checks 

The correlators Ily^ given in Eqns. (^, ^4|) involve the finite functions Afin, B defined in 



the Appendix as well as the renormalized coupling constants Lf^ ^Q, Hi'\Ov, Pv,Qvj Rv- 
Note that upon taking the (1+0) component, Oy and the direct dependence on L^^ drop 
out (the indirect, through n, remains). Our result satisfies the following consistency checks: 

• The low energy constants and the function Afin (M|,) depend on the renormalization 
scale n- Upon using Eqns. (|25| - p2D as well as the renormalization group equations for 



the pertinent low energy constants it can be shown that the sum of all contributions 
is scale independent. 

The SU(3) limit of equal quark masses provides another check: The J = compo- 
nent vanishes identically, as it is in proportion to the difference of two meson masses 
squared. The J = 1 component reduces to the known results for the isospin or hy- 
percharge component in the degenerate case. Although non-trivial, this consistency 
check concerns the SU(3) conserving part of the amplitude only. 



Going away from SU{3), the AdemoUo Gatto theorem is satisfied: The SU{3) 
breaking effect appears only in the second order in the quark masses, i.e. oc (m^ — m)^. 

Individual spin components J = 1, exhibit poles at = 0. These poles must be of 
kinematical origin since there is no single particle state present which could generate 
such a singularity. Using - instead of (^OD - the new decomposition 
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(0)/ 



(41) 



both new components, Ily^'^\q^) and q^Ily\{q^) are found to be regular at = 0. It is 
thus more convenient to work with these combinations, in particular when considering 
dispersive representations based on analyticity assumptions. 



The explicit form (p3| , of the vector correlator allows us to extract the corresponding 
spectral functions pyl = ^Imlly] for J = 0, 1 in a straight forward manner: The 
absorbtive parts are contained exclusively in the functions -B(g^,M^,M|), i.e. 



1 iy{q'^,Mf,Ml 



IGir q^ 

2 ^T\J , T\J \2 



^'eiq'-iM^ + M^y) 



with 



iy{q\ Mi, M^) = q' - {M, + 



1/2 



g2 - (Ml - M^y 



1/2 



(42) 
(43) 



By simply rearranging contributions one observes that the result for J 
expressed in terms of the on-shell vector form factor ([T5|), viz: 



1 can be 



Pv,s 



647r2 



u\q^Ml,M'^: 
q6 



e{q' - (Mk + M^Y) + {71 ^ri) 



Likewise, but with considerable more effort, it can be shown that the J = spectral 
function is given in terms of the on-shell scalar form factors (|r^, ^), viz: 



(0) _ 3(Mi - M2 

Pv,s — ' 



2\2 



v{q 



2 Ml, Ml) 



e{q' - {Mk + M^f) . \ft{q')? + (vr - r/) . (45) 



647r2 



In deriving 



we have used the GMO relation as well as the expansion 



|/;;o«r = l + 2ReA/;;o^ + 0(p' 



px 



(46) 



which follows from f\.'^ = 1 + ReA/+;5^ + ilmA Eqns. (|4|, H) are of course 
nothing else than the unitarity conditions for the absorbtive parts to the order we 
are interested in. Phenomenologically, the XPT expressions for the J = 1 spectral 
function may (at best) be used in the immediate threshold region only, since the 
prominent resonance K*{%92) is very close by. 

The physical content of our result can be further assessed by studying the asymptotic 
behavior ^ oo by means of the operator product expansion. Accordingly, the dispersion 



relations for Iiy'].'^\q^) and q^Iiy'g{q'^) need at least one subtraction. In the chiral represen- 
tation this fact is reflected by the occurrence of contact terms, e.g. those in proportion to 



r(o), 



nf'' and C^i . They depend on the particular way the ultraviolet divergences have been 
regularized. These contact terms are seen to disappear from our result upon either taking 
derivatives with respect to q^ or, in the case of the J = 1 + combination, by forming 
flavor breaking differences of the flavor components nyj°^(g^) with / G {3, 8, s}. Hence, the 
proper physical observables are 



.(0) 
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<f (g^), TtWM^K^')' ^=1,2,... (47) 



and <+°)(g2)-ni^j,o)(g2)^ A, G {3, 8, . (48) 

These functions depend on only two counterterm coupling constants of the 0{p^) chiral 
Lagrangian, i.e. PvifJ') and QvifJ')- Since the constant PvifJ') has already been determined 
to high accuracy |jl3[, we concentrate in the following on the flavor breaking differences (^8]). 



B. A low energy theorem 

The general structure of the flavor breaking differences is 

4;j:°)(g^)^nl^;;)(g^)-nS;j°)(g^) 

= D};)ri^") + Dllr-'^^'ig^ /i) + DJ^%{q\ /i) . (49) 

The first term in the last line of (^) denotes the total one-loop contribution. The second 
and third term represent the two-loop contribution which is separated into finite loop- and 
the finite local 0{p^) counterterm pieces, respectively. Explicit expressions for these terms 
can be obtained from (|33| , and Ref. Note that the 0{p'^) counterterms L^q, h[^^ 



drop out in such differences. Hence, the only terms containing free parameters are 



(Mi-Mf 







16 (Ml - M^' 
F2 



Consequently, there is one combination of vector-current correlators which is free of any 
counterterm coupling constant, viz. 

AUy{q') ^ <r)(g^) + 3<r)(g2) - 4<r)(g^) = AU\,-'^^' {q') + Atf^-^°°P(g2) . (51) 



Eqn. (|5T| ) is a low energy theorem valid at the two-loop level. It shows that despite the 
many low energy constants occurring at this order of the low energy expansion it is possible 
to give parameter free predictions for suitably chosen combinations of observables. 
The RHS of (|51|) is readily evaluated numerically. At = 0, e.g., we obtain 

Anv(O) = 0.00513 + 0.00125 , (52) 

where the 1-loop and 2-loop contributions have been displayed separately. We observe a 
correction relative to the leading order term of ~ 25%. This prediction for Any(O) can be 
tested if the vector spectral functions Py/^"* are known to high accuracy. For instance, the 
unsubtracted dispersion relation for AUyiq'^) implies the chiral sum rule 

Ally{0)= ds^_Y:^^lPYl (53) 

J4M2 S 
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Saturating the spectral functions with narrow resonances for p(770), 07(782), 0(1020) and 
K*{892) gives the estimate 

Any(O) ^ 0.039 + 0.085 - 0.121 = 0.003 . (54) 

The high degree of cancellation present in (|5^) clearly calls for the inclusion of finite width 
effects, higher resonances as well as multiparticle contributions to the spectral strengths. 
However, currently available data for hypercharge and strange components of the vector 
spectral functions do not prove accurate enough to make this test conclusive. We shall 
present a more detailed analysis involving the difference Pvt^^ ~ Pvt^^ in section V below. 



V. INVERSE MOMENT FINITE ENERGY SUM RULES AND 
DETERMINATION OF Qy 

The two-loop representation of Hyf^ given in the previous sections can be combined 
with the isospin or hypercharge component of the vector current correlators to form 

physically observable quantities. In Ref. [0 the difference between correlators of isospin 
and hypercharge currents was employed to estimate the finite counterterm Qv- It has been 
pointed out later that isospin violating corrections have a rather large effect on the outcome 
of that analysis fl^ . The reason is that individual components get mostly canceled such that 



the isospin breaking effect on the difference is enhanced. Isospin violation affects primarily 
the hypercharge component and is hard to be treated in a model-independent way. The 
aim of this section is to present an analysis based on the difference of isospin and strange 
components of the vector current which is less sensitive to such corrections. Thus, we shall 
consider the flavor breaking difference 

n(g2)^ny,r)(g2)-nS;j°)(g^), (55) 

which is free of kinematical singularities at the origin. The physical content of the XPT 
(chiral perturbation theory) representation of such combinations is conveniently analyzed 
by means of inverse moment finite energy sum rules (IMFESR) [l^Jl^. Unitarity and 
analyticity imply that the correlator ([55|) satisfies 



^ds^J^ = 0,n = OA,... (56) 

where C is the contour shown in Fig. | and w{s) is a conveniently chosen weight function 
which is analytic inside this contour. In the following we shall consider polynomial weight 
functions with w{0) = 1. The contribution from the small circle can be computed using the 
XPT representation of the correlator, i.e. 

1 ,^ , If , n(s)w(s) f^o p(s)w(s) ^ ^ . 

(nw;)(0) = — ds + / ds I J , n = 0,l,... (57) 



n!(cig2)nv 27ri y|s|=so s"+i 5"+^ 

with p{s) being the spectral function of the correlator under consideration, i.e. 

p{s) = -Imni^f ) is)e{s - AMI) - -Im4^^°) {s)e{s - {Mk + M^f) . (58) 
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The first term on the right hand side of (|57|) can reliably be evaluated using the operator 
product expansion for the correlator, provided the radius of the circle, Sq, is chosen large 
enough. By contrast, the second term on the RHS of { ^7\ ) needs the difference of hadronic 
isospin and strange vector spectral functions as input. The only source of such hadronic 
spectral functions is provided by the tau decay data [^|11|; for this reason sq cannot be 
taken much larger than the tau mass squared. 

The weight function w{s) needs to satisfy 3 more conditions to make the calculation 
feasible. First, w{s) must be chosen in such a way that the perturbative series in of the 
leading OPE contributions is well controlled. Second, the region close to the real axis, s = Sq, 
must be suppressed since the OPE is not reliable in that region. Weight functions without 
this suppression substantially violate local duality. In practice, this suppression is highly 
welcome for reasons which relate to the third condition: As s approaches the tau mass, 
data become more and more uncertain. As a consequence, the integral over the hadronic 
spectral function is afflicted with large error bars — unless the weight function suppresses 
the contributions from, say, s > (1.4GeV)^. It is this point which makes the method of 



IMFESR superior to the conventional inverse moment chiral sum rules |T3| , pT| -|23|] : The high 
energy behavior of the weight function can be worse than what a once subtracted dispersion 
relation would require. Nevertheless, as long as the weight function exhibits one or two zeros 
at s = So, local duality is extremely well satisfied, even at scales below s = [Q. 



In the following, we shall study the IMFESR (pTQ for n = using the weight-function 

1 s 

w{s) = {l-xf-{l+x+-x'^), x = — . (59) 

2 So 

This weight-function satisfies the conditions listed above. The factor (1 — x)^ proves efficient 
in cutting off the region above 2GeV^. 

Next, we specify the input needed to evaluate the inverse moment finite energy sum rule. 
We shall consider in turn the XPT, the OPE, and the hadronic side of the IMFESR. 



A. Chiral constraints at = 

The correlator at = is obtained by straight forward evaluation of the representation 
of Ily-^^-* given in section III of this paper and Ily-g"'' as calculated in Ref. 0. Some care 
has to be taken in order to show the cancellation of singular terms in the limit — > 0. 
Employing input parameters = 140 MeV, Mk = 495 MeV Fq = 93 MeV, and, at 
renormalization scale fx = rup, = 0.0014 ± 0.0005 [|5|, = 0.00678 ± 0.00015, 



= -0.00513 ± 0.00019 [|T9l we obtain numerically 



4 



n(0) = 0.0053 + 0.0014 + -=;^{M'k - M^)Qv{m;) . (60) 



The first and second term on the RHS of (0) are finite one- and two-loop terms, respectively. 
The counterterm coupling constant Qy has been estimated previously, with the result 



Qvifi = rUp) = (3.7 ± 2.0) 10"^ [13] 
Qvifi = mp) = (3.3 ± 0.4) 10'^ 



(61) 
(62) 
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where the second number has been obtained with a phenomenological ansatz for the 38- 
component of the vector spectral function. Using the first of these values, the last term in 
(|60|) yields numerically 0.0038 ± 0.0021, i.e. at scale /i = rup it is the dominant effect of the 
two-loop contributions. 



B. OPE of the correlator 



The contributions of dimension D = 2,4, 6 are known from previous work [pq-|28[|. The 
D=2 operators are just the flavor breaking quark mass corrections: 



2\D=2 



n(Q 



where 



1 3 
02 4^^ 



, 7 
1 H — a, 
3 " 



+ 19.9332a: + 0(< 



-as + Oia^ 



+ 0{ml] 



(63) 



-q^ denotes the Euclidean momentum, and = as(<5^)/7r. Since terms in 



proportion to msniu are further suppressed by a^, 
excellent approximation. 

m 



keeping the mj. term only provides an 



Following the prescription of Ref. 
contour with radius sq 



the integration of this function over the circular 
is performed numerically. We use the 4-loop (3- and 7-functions 
|, |3T|] for the running of mass and strong coupling constant. Weighting the correlator with 
w{s) given in (^91), we thus obtain 



B 



D=2 



1 

2"7ri 



/ ds 

s\=so 



U(s)w(s) 



47r2 



mt 



[(1.634 + 0.303 + 0.177) ±0.088] . (64) 



The three terms in parenthesis represent the contributions at 0(1), 0(0^), O(a^) in the con- 
tour improved perturbation series. The error reflects the uncertainty due to terms of order 



assuming a geometric growth of the series in at Q 



The particular combina- 



tion of moments appearing in the weight function w{s) makes sure that even the resummed 
geometric series yields corrections not larger than the size of the term. 

In evaluating Eqn. (^) the four-loop beta function is used, but values of extracted are 
usually somewhat dependent on which order of perturbation theory is used. To get a handle 
on the uncertainty brought by the truncation in the perturbative series, we mention that 
the bracket in Eqn. (|6^ ) is changed into [1.668 + 0.316 + 0.190 ± 0.096] if two-loop running 
is used. In other words: [2.174 ± 0.096] instead of [2.114 ± 0.088] is a 3% effect. 

As a result, the total error is dominated by the uncertainty of present evaluations of 
the strange quark mass. For a numerical evaluation we use m<j(/i = m,-) = (130 ± 25)MeV 
2§,|23,|g, leading to 5°=^ = (8.5 ± 3.3) 10^1 This corresponds to ^ 16% of the leading 
order one- loop term in (|60D. 

The D=4 operators consist of two types, the quark mass corrections and terms in pro- 



portion to the quark condensates. The leading quark mass correction is in proportion to 
a~^m% Despite th 
is of the order 10" 
above. The condensate terms read 



5 ^ enhancement, the net contribution to the RHS of the sum rule 
i.e. completely negligible compared to the D=2 contributions found 



n(g 



2^D=4 

condensate 



-Ims , ^ _ , 
-TTT^—imuu) 



59 13 



(65) 
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The ratio Tc = {ss) / {uu) lies in the range 0.7 < < 1 
we employ leading order XPT expressions, i.e. nis/rfi = 
this manner, the contribution to the IMFESR under consideration is also found to be of the 
order 10~®. Hence, we can safely neglect the contributions of D > 4 operators in the OPE. 



2^j27[] . For numerical evaluation 
25.9 and {muu) = -FlMl/2. In 



C. Hadronic spectral function 

The isovector- and strange component of the vector spectral functions are constructed 
as follows: i) close to threshold we employ the XPT two-loop expressions given in and 
this paper. These results depend on the low energy constants L^,Lg. However, due to the 
smallness of the spectral functions in the threshold region compared to the resonance region, 
uncertainties in these LECs have no detectable effect on the IMFESR (0) with n=0. ii) 
above the threshold region and up to the tau mass we employ the data as presented by the 
ALEPH collaboration [^,|n|]. We discuss this input as well as the corresponding evaluation 
of integrals over hadronic spectral functions for isovector- and strange component in turn: 



isovector vector spectral function: the impressive set of ALEPH data can be used 
over the whole range of integration. Data are accurate to ~ 3 % for s < 2GeV^. 
The region above this scale is however strongly suppressed when integrating with the 
weight function (^9|). Table 1 shows the hadronic integrals 

Byt (Smax, So)= / ds — (66) 



for various components / of the vector spectral function and different upper integration 
ends Smax- The parameter sq refers to the weight function (^) and is taken as large 
as possible, i.e. sq = ml ~ (1.777 GeV)^. The first line in Table 1 is split into three 
columns to illustrate the effect of not extending the integration range in (^) to its 
upper end sq, but to 1 or 2 GeV^ only (while keeping the weight function ( pUD with its 
zero at sq). In the isovector case, the dominant contribution to the hadronic integral 
stems from the region below ~ 1 GeV^, where the spectral function is still dominated 
by individual resonances. The error in is predominantly due to the correlated 
part of the error in the individual bins, i.e. mostly due to the uncertainty of the r 
branching ratios. We are working in the isospin limit where the J = component of 
the isovector vector spectral function vanishes. 

(b) strange component of the vector spectral function, J = 1 + 0: in a recent publication 
[pj]] the ALEPH collaboration has presented results for the strange component of 
vector plus axial-vector spectral functions obtained from tau decay data. Unlike in 
the isovector full separation of vector and axial-vector components was not 

achieved. However, as we shall argue below, the information is sufficient to evaluate 
the difference of hadronic integrals we are seeking to an accuracy of ~ 17 %. To this 
end we analyze the data according to the hadronic state X in the decay t ^ + X: 

— X = Ktt: the decay to the Kn final state is a pure vector current transition. The 
main feature of this component is the K*{892) resonance, but there is a minor 
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contribution associated with the interference of i^*(1410). The resulting integral 
is given in the third line of Table 1. Again, the error is primarily due to the 
uncertainty in the overall normalization. 

X = K2ti: the final state with two pions is due to both, vector- and axial- 
vector transitions. As detailed in []TT| the corresponding spectral function can 



be understood in terms of the two axial- vector resonances i^'i(1270), i^i(1400) 
and the vector i^*(1410). The mass resolution and statistics are not sufficient 



to separate the i^i(1400) and K*(1A1Q). We follow [jTT| and use an effective 
resonance instead, averaging the parameters of the two states. The vector part 
is then estimated by fitting p^^] + p^^] to two resonances with masses mi = 1.270 
GeV, 7712 = 1-405 GeV. Switching off the axial ii'i(1270) yields an estimate for 
the contribution from the resonances at s = (1.4 GeV)^. Half of the remaining 
spectral strength is assumed to be vector, with an estimated error of ~ 50 %. 
This procedure is roughly consistent with the vector- axial-vector separation of 
branching ratios shown in Table 7 of Ref . . The contribution to the hadronic 
integral B^2tt is roughly 7 % of the Ktc component, see Table 1. 

— X = Krj and K3tt: this component has not been separated into vector- and 
axial-vector contribution. However, the spectral strength associated with these 
final states is relevant only above s ~ (1.3 GeV)^. Due to the suppression of 
this energy region by the weight function, the net contribution to the hadronic 
integral is small. We assign half of this component to the vector current and 
employ the rest of the spectral strength as an estimate for the error. 

— X = KniT for n > 4: these components are small and non-zero only above 
s = 2 GeV^. For this reason they can safely be neglected. 

Again, the different columns in Table 1 illustrate the effect of not extending the in- 
tegration range in (36) to its upper end Sq, but to 1 or 2 GeV^ only (while keeping 



the weight function ( |55| ) with its zero at sq). The contributions to the strange vector 
spectral function are not dominated by the low-energy part as impressively as the 
non-strange part, but still the integrals are almost saturated at ~ 1 GeV^. 

strange component of the vector spectral function, J = 0: the main contribution to 
the exclusive J = component is provided by production of the scalar resonance 
i^'o(1430), which decays almost exclusively to the Ktt final state |TT[]. In order to 



obtain an estimate for this component we use an Omnes representation as given in 
Ref. 1 5^. Employing also input parameters as detailed in yields the integrals as 



shown in Table 1. The error is an educated guess of the various uncertainties entering 
the Omnes representation. This effect is already accounted for in the Kir component 
and should not be included when summing the strange contributions to By^^. 

The hadronic integrals summed over the strange components of the vector spectral func- 
tions is given in the third to last line of Table 1. Subtracting this result from the hadronic 
integral over the isovector component and adding errors in quadrature we finally obtain 

B^f - B^"^^ = 0.0087 ± 0.0013 . (67) 
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Individual flavor components thus have canceled to ~ 60 %. The error stems mostly from 
the error on the /TTT-component of the strangeness carrying vector spectral function, but the 
error on the ii'27r-component is non-negligible as well. 

The final step is plugging the results of Eqns. 0, into the IMFESR (|57D for 
n = 0, which yields 

Qy(/i = mp) = (2.8 ± 1.3) • 10"^ (68) 

This result is consistent with the determinations ( pT| , |B2|), but with the advantage of being 
both, relatively precise and free of model assumptions. 

The variation of Qv with renormalization scale fi is known 0: 



Qvi^^) - Qv(/xo) = (Lf\M,) + Lffi(M,)) In ^ . (69) 



Qv is thus strongly varying with scale. Explicitly we obtain e.g. Qy(500MeV) = (4.2 ± 
1.4)10"^ and Qvi^GeV) = (2.0 ± 1.4)10"^. However, the method employed here deals with 
RG-invariant quantities, and the variation with scale does not lead to any further ambiguities 
in the determination of Qv- We also note that the coupling Qv is the same in both, the 
A-subtraction scheme used here and the MS-scheme [jl5| . 



In our determination, sq was chosen as large as possible by having it equal to m^. Of 
course one may ask whether the total result is stable under a variation of the radius of the 
outer circle in Fig. ^. Since there are no data from r-decays with s > m^, the only alternative 
is a radius which is smaller than our choice sq = m"^. Concerning the contribution from the 
hadronic spectral integral, there is a clear numerical answer to this question. Repeating 
our calculation, now replacing [both, in the upper boundary of the integral (^) and in 
the weight function (|5|)] sq = (1.777 GeV)^ by (1.555 GeV)^ and (1.333 GeV)^, we find 
the values collected in Table 2. As one can see, the total hadronic spectral integral (3 — s 
component) is extremely stable under this shift. On the other hand any numerical statement 
about the shift in the contribution from the outer circle is somewhat questionable. Lowering 
the radius in the integral ( |64D to |s| = Si = (1.555 GeV)^ results in the bracket in the r.h.s. 



being replaced by [1.570 + 0.302 + 0.180 ±0.090]. Going further down to a radius as small as 
|s| = S2 = (1.333 GeV)2 gives a bracket which takes the value [1.474 + 0.296 + 0.181 + 0.093]. 
The corresponding shift in B^^"^ is small compared to the error brought by the uncertainty in 
present evaluations of the strange quark mass. Hence, the contribution from the outer circle 
seems to be insensitive to the precise value of the radius too. The real problem, however, is 
that the expansion in contributions at different orders in the OPE breaks down at s-values 
substantially below 4GeV^. In this respect choosing sq = amounts to stretching the 
validity of the perturbative expansion to the very end. All one can say is that our result is 
stable under small variations of the outer radius in Fig. |. 

Regarding the error, it is clear that a more accurate determination of both the strange 
and the non-strange r branching ratios as well as a more accurate measurement of the Kn- 
component of the vector spectral function below s = (1.4 GeV)^, in particular in the i^'*(892) 
region, will transform into a more precise determination Qv- As for the K27r component, a 
full separation into V- and A- part would further diminish the error. 



Our final value for Qv supports the claim by Maltman and Wolfe [0 that isospin 
breaking effects tend to reduce the estimate (|6lD. 
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VI. CONCLUSIONS 



We have calculated the two-point function of strangeness carrying vector currents to 
two- loop standard chiral perturbation theory, in the isospin limit. Our fully renormalized 
finite result passes several consistency checks: The sum of all contributions does not depend 
on the renormahzation scale In the SU{3) limit of equal quark masses it reduces to the 
known resuhs for the two-point functions of isospin and hypercharge currents. The SU{3) 
breaking part satisfies the constraint imposed by the AdemoUo Gatto theorem. Further- 
more, kinematical singularities at zero momentum transfer present in the spin and spin 1 
components are shown to be absent if a convenient and well-known decomposition is cho- 
sen. Last but not least, the imaginary parts of the correlators for spin and 1 obey the 
unitarity-conditions along the two-particle cut implying that they can be expressed in terms 
of the vector and scalar meson form-factors. 

The calculated correlators 11^] for J=0,1 depend on several low energy constants (LEG). 
Two of these constants are related to contact terms. In physical observables, which do not 
depend on contact terms, only two low energy constants of the 0(p^) chiral Lagrangian 
remain, i.e. Py and Qy- We have shown that at 0{p^) in the low energy expansion the 
combination (llyjl^'^'' -|- Slly^^'^'* — 4ny^^'^'')(g^) is free of any unknown LEG and consequently 
predicted in terms of pseudoscalar meson masses and decay constants only. This low energy 
theorem (LET) can be tested once the vector spectral functions pv,3, Pv,8 and pv,s are known 
to high accuracy. Presently, data are not accurate enough to make such a test conclusive, 
but the situation might change within a couple of years. 

Finally we have presented a determination of the low energy constant Qy based on a 
inverse moment finite energy sum rule. The sum rule employs the difference of isospin and 
strange components of vector spectral functions and is therefore expected to be less sensitive 
to isospin breaking effects than previous attempts. Our analysis is based on the recently 
released ALEPH data of r-decays into strangeness carrying final states. Presently, the error 
bars on the data are still too large to make the uncertainty of our first-principle based Qv 
determination smaller than the statistical error of previous model-based calculations. How- 
ever, the fact that the main uncertainty stems from the dominant decay channel means that 
our hope for more precise experimental data might be realistic. The very same experimental 
effort would also be welcome in view of the low energy theorem discussed above. 
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APPENDIX: SURVEY OF ONE-LOOP INTEGRALS 



Here we give a survey of the integrals occurring in the computation of the formfactor and 
the desired correlator. We shall keep d arbitrary and isolate the divergences near (i — >■ 4. 
The simplest quantity is the scalar integral 

A{m^) = f — i— (70) 

which, by standard means, is evaluated to give 

^('"') = ^r(i-5)K)""'-". (71) 

From this expression it is found to split 

A{m'^) = -2irn^ fi^'^-^^ A + Afin(m=^) with (72) 



where ( [74D is the part of A which stays finite under d ^ 4 in this so-called A-scheme. 
Each of the contributions in (|72D depends on the artificial scale fi and the prefactor /x^'^"^) is 
introduced in order to give them the appropriate mass dimension in d space-time dimensions. 
We shall define 



dH 


1 






{2'nY (^^ - ml) 


((£- 


g)2 


- ml) 


dH 


% 






(27r)^ (£2 - ml) 


((£- 


g)2 


- ml) 


dH 









' = y (2^ {P- ml) {{I-,)- -ml) ^^'^ 
B^{q,ml,ml) = [ -^J-— (76) 



m„ = y ^ ii^-ml)iii-q)^-ml) ^^^^ 
as well as the finite quantity 

B{q\ ml ml) = B{ql ml ml) - 5(0, ml ml) (78) 
which obviously satisfies 5(0, m^, 7712) = 0. Then the exact relations (in d dimensions) 



,2 ^2^ ^ (^-2) A(m 2 
2 m2 
2 2, A{ml) - A{ml) 



5(0,m^m2)= (79) 



5(0, ml ml) = ^ (80) 



mf — m: 



follow by taking a derivative in ([70| , |7ID and by partial fraction decomposition in ( [TSl 
Invoking covariance, 5^ and 5^i, in Eqns. (|7^, |7^) may be decomposed as 
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BM.mlrnl) = q,B,{q\mlml) (81) 
B,M.mlml) = q,q^B,,{q\mlml)+g,^B,,{q\mlml) . (82) 



Considering g^-B^ one finds 

B^{q\ ml ml) = ^ (1 + B{q\ m^ m^) - p^^" (83) 



1 ml- ml 2 2 (A(mf) - A {ml)) 
2(1 + — — ) , "^1, ^2) ^ 



-[l + ^—^)B{q\mlml) + ^^-^, \^ 84 



where in ( p4D one would have to use Eqns. ([79| , |80|) for performing the equal mass limit. 
Considering q^B^y and g^vB^y one finds 

(<i - 1) B,. = ^ (1 + ^^i^) B, - !^ B + ±^ ^(m^) (85) 

ci(g^ + m\ — mlY — Aq^m\ — 



4g2 

{^-d)m\-d{q^ -ml) ^ d{q^ - mj) - {A - d)ml . 

4(mf - ml)q^ ^^"^'^ 4(m? - ml)q^ ^^""^^ ^^^^ 

(c/ - 1) = -]^{q^ + m\-ml)Bi+ m\B-]^ A{ml) (87) 

q'^ — 2q^{m\ + m^) + (mf — 7/12)^ — 

g2 _ 3^2 _ ^2 g2 _ ^2 _ 3^2 

4(mf-mi) ^^"^'^^ A{ml-ml) "^^"^'^ 

where in (^6]) and (|88D one would have to use Eqns. (|79| , pOD for performing the equal mass 
limit. A combination frequently occurring in loop integrals is 

T^,(g, ml ml) = AB^, - 2q^By - 2g,5^ + q^q.B (89) 
= 9^.AB22 + qM^B^i - ABi + B) (90) 
= PfJ ■ {aB,2 + (4^21 - 45i + B) q') + P« ■ (4522) (91) 

where, in the last line, we have introduced the projection operators 

P^'J = '-f and pa) = ,,,-M!l. (92) 

Denoting T^o) = 4^22 + (4^21 - 45i + B) q^ and T^^) = 4^22 the spin-0 and spin-1 compo- 
nents, respectively, and using the decompositions (Q ^ ^Sj) along with ( [78| , [79| , ^OD these 
two pieces are found to split up into divergent and finite parts 

tIS (ql rnl ml) = -2i (m? + ml) A (93) 

(g^ m?, ml) = rnl) + AU<) + A,^{ml) (94) 

iq', rnl rnl) = | (g^ - 3(m? + m^)) /i^^-^) A (95) 
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6 ^ q'^ 



q"^ — 3m? — mi, , ,9, q"^ — mi — Smi , , n, 

+ 2 2 ^fin "^I - 2^ ^ ^fin 

mf — m2 mf — m2 

^(g2-3K + m^))}, (96) 
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where terms contributing to 0{d — 4) have been dropped and where the equal-mass hmit is 
trivial for every quantity but Tj^]^ , for which it takes the form 



T£)(g^m^m^) = -i{(g^ -4m^) S(g^m^m^) + {q'-6m')^^^ - ^} . (97) 



TABLES 



'5max(^ "^o) • 


IGeV^ 2GeV^ 


So = (1.777 GeV)2 


3-component, J=(l) = 


(1+0) 0.0240 ± 0.0005 0.0254 ± 0.0005 


0.0257 + 0.0005 


i^vr-component 


0.0146 ± 0.0009 0.0155 ± 0.0010 


0.0155 + 0.0010 


K2tt- comp onent 


0.0000 ± 0.0002 0.0008 ± 0.0004 


0.0011 + 0.0006 


Kt]- and i\:37r-component 0.0000 ± 0.0002 0.0003 ± 0.0003 


0.0004 + 0.0004 


sum s-component, J= 


(1+0) 0.0146 ± 0.0009 0.0166 ± 0.0011 


0.0170 + 0.0012 


3-s, J=(l+0) 




0.0087 + 0.0013 


s-component, J=(0) 




0.0011 + 0.0004 



Table I: Hadronic integrals -Bvf/'"^(smax, sq) as deGned in for various components f of 
the total vector spectral function py^ as defined in Eqn. ( ^ ). Weight function is (|5^ ) with 
SQ = ml = (1.777 GeV)2. 







W = S2 = (1.333GeV)^ 


Smax = si = (1.555GeV)^ 


3-component, J=(l) = 


(1+0) 


0.0160 + 0.0004 


0.0213 + 0.0005 


i^'TT-component 




0.0073 + 0.0005 


0.0119 + 0.0008 


ii'27r-component 




0.0000 + 0.0001 


0.0003 + 0.0003 


Kt]- and ir37r-component 


0.0000 + 0.0001 


0.0001 + 0.0002 


sum s-component, J=l 


:i+o) 


0.0073 + 0.0005 


0.0123 + 0.0009 


3-s, J=(l+0) 




0.0087 + 0.0006 


0.0090 + 0.0010 


s-component, J=(0) 




0.0003 + 0.0001 


0.0007 + 0.0003 



Table II: Hadronic integrals B^j''^ {s2, S2) as defined in (|6^ ) for various components f of 
the total vector spectral function pyj as defined in Eqn. ( ^) . Weight function is (|5^) with 
So si = (1.555 GeV)^ and sq ^ S2 = (1.333 GeV)^, respectively. Either one of the two 
columns corresponds to the rightmost column in Table I. 
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FIGURES 




(a) 



(b) 



(c) 



FIG. 1. One-loop Fcynman graphs for Ilv,s- Dots and crossed circles denote vertices from the 
lagrangian of order and p'^ respectively. Wavy lines are external vector currents, straight lines 
denote Goldstone boson propagators. 




(a) (b) (c) (d) 





(e) 



(f) 



(g) 





(h) 



(i) 



a) 



FIG. 2. Two-loop Feynman graphs for Ilv,s- 1^ addition to the conventions as in Fig. 1, the 
crossed square denotes a vertex from the lagrangian of order p^. Diagrams (g) and (h) shall 
comprise the corresponding modification to the incoming vector-vertex as well. 
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(a) (b) (c) 




(d) (e) (f) 

FIG. 3. Feynman graphs for meson form factors at one-loop. The partial wavefunction renor- 
malization convention used in the calculation affects the weight of diagrams (b, c) - see text. 





Im(s) 




c / 














— t , V 


1 F 









(s) 



FIG. 4. Integration contour for the inverse moment Rnite energy sum rule (|5^). The cut on the 
real axis starts at Sth and extends to infinity. The dot at the origin indicates possible poles in the 
case of inverse moment finite energy sum rules. 
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